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ABSTRACT 

We present simple inequalities for the Riemann problem for a Hamilton-Jacobi 
equation in N space dimension when neither the initial data nor the Hamiltonian 
need be convex (or concave). The initial data is globally continuous, affine in each 
orthant, with a possible jump in normal derivative across each coordinate plane, 
Xi = 0. The inequalities become equalities wherever a “maxmin” equals a “minmax” 
and thus an exact closed form solution to this problem is then obtained. 
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We shall be concerned with solutions to the following differential equation 
(H-J) (ft + H{D x fi) = 0 in R N x (0, oo) 

where HsC(R N ) and D x <fi = {fix i , . . • , fi Xn ) is the spatial gradient of fi. 

We shall take special (Riemann) initial data. Let uf,u~ be constants and 


define 

tii(x) = uf if X{ > 0 

(1) 


u,(x) = u~ if n < 0 

for i = l,. 

. . , N. Then take: 


N 

(2) 

fio(x) = A + u i( x ) = A + 


»=i 


We wish to solve H-J with initial data (2) in the class of viscosity solutions 
as defined in [2]. The four properties of viscosity solutions that we shall need here 
(apart from existence, uniqueness derived in [2]) are: 

(Pi) The solution fi(x,t) is a non-decreasing function of the initial data. 

(P2) The partial derivatives < p x . satisfy a maximum principle at points of continuity, 
i.e. for * = 1, ... , N: 

min(u~ ,uf) < <fi x . < max(u~, «t). 

(P3) The speed of propagation is finite. 

(P4) If ifi(x 2 , . . . is a viscosity solution of 

4>t 4 = 0 
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for a constant t>i then 


V>{x,t) = V\X\ + ) 

is a viscosity solution to (H-J). 

Let 

Q = fli x x • • * x &N 

where for i = 1, . . . , N: 

fli = { v / win(u~,uf) < v < max(u~, it+)} 

Finally we let, for i = 1, . . . , N: 

Xi = sign ( uf - u~) 


For convenience only, we order the indices so that 

Xi = i = 1,2, . . . ,j 
Xi = ~ 1, i = j + l, . . . ,N 

(j might be 0 or N + 1). 

We now state: 


THEOREM 1. The viscosity solution to (H-J) with initial data (2) satisfies : 
A + max max • • • max min • • • min [x • v — tH(v )] 

i>i*Oi 172*02 VjeQj v ; -+i*0j+i ujv*0iv 

(3) ^ ?(*» *) ^ 

A + min ••• min max • ■ • max[x • v — tff(v)] 


We note that if all the Xi — then this solution is just max[x • v — tH{v)\, if all 

17*0 

the Xi = —1 it is min[x • v — tff(v)]. Otherwise we have a pointwise inequality which 

v*0 
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gives the exact solution whenever the first and last terms in (3) are equal. This 
occurs e.g. if H{y ) = Hi(v \, . . . , Vj) + HiiVj+i * . . . , i*e., if the Hamiltonian 

separates, and in many other cases. 

The rest of this paper consists of the proof of this theorem, and some remarks 
about both conservation laws and numerical approximations to H-J. 

It is easy to see that the solution to the Cauchy problem satisfies 

(4) <p(x,t) = tg(j) + A = tg(Q + A 
where g satisfies: 

(5) g = < • D ( g - H(D ( g ) = ‘(sifO. j 2 (0, • • ■ , Jw(0) 

where D^g is continuous. 

In H-J, we let r = t, y f - = x, — £2 for £ fixed. H-J becomes 

v? r + H(Dy<j>) - CD 9 <p 

(H-J 1 ) 

= <p r + H l (D y y) — 0 

with the same initial data (2). 


Thus, by (5), to evaluate g(() we need only evaluate —H l {D y g) at y = 0, for 
any t > 0. From (Pi) above we know that (D y g) y= o lies in 12 for t > 0. Moreover, 
if we integrate (H-J) 1 from t = 0 to t = At we have 
V?(0,At) = A- AtH l ((D y g) y -o ) 

— \DF<p o(0),I>r V o(0)) 

Here: 


(7) 


DSrtO ). ± (*(**«) 
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where c% = {0,0, ... ,1,0, the i 1 * 1 unit vector, and ; u 2 , u 2 , • • • , u pj , 

is determined by (6). 

This formula can be interpreted as a numerical algorithm. Suppose we are 
given a grid 

x )i ~ • • • » -W; Ji = ■ ■ • 

and values of a discrete function xf)j = V’iiia—iiv Then for each j , we construct the 
piecewise affine function which, in each of the 2^ orthants centered at j , interpolates 
ipj and its N nearest neighbors, i pj± €i for t = 1, . . . ,N. From (P3) above, if 

( 8) (cfl) 

iml N 

where is the same as ft with each u~, uf replaced by D X J rpj, £>+>;, then the 
solution to the initial value problem (H-J) 1 with the above affine initial data in the 
diamond centered at j when evaluated at x — Xj and t = At is independent of the 
values of the initial data outside of this diamond. 

Thus (6) (with y? o (0) replaced by ipf and y>(0, At) by V>" +1 ), gives us a mono- 
tone finite difference scheme approximating (H-J) 1 which is in differenced form with 
numerical Hamiltonian H 1 . These concepts were introduced in [3]. The scheme is 
monotone, which means that the right side of (6) is an increasing function of all the 
because of property (PI). The function H 1 is called Godunov’s Hamiltonian 
by analogy with the definition of Godunov’s scheme for conservation laws in one 
space dimension [5]. The scheme is consistent, which means 

H 1 (u i ,ui; u 2 , 112 ; • • • ; U N ) = H x («i, U 2 , . • • , «n) 


Monotonicity implies that 

H (u^" , Uj J U 2 , U 2 ,,..., Ujy) 
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is a nonincreasing function of all the uf and a non-decreasing function of all the 
uf . In particular, for N = 1, this means for any vjefi = j : 


(9) sgn (uf - - tfVi)] 

= sgn (uf - v x )[H X (u^ , -uf) - H 1 (v! , uf )] 

+ sgn Oi - uf )[fi' 1 (ui,uf ) - tf 1 ^!,^)] 

< 0 

But, by (P2), H l (u+ ,u^) = ff 1 (u 1 ) for some ui in fl. Thus we have 

(10) £>(«+ uf) = X i minxi^i) 


(This formula was obtained earlier in [6]). Now we proceed inductively. Suppose, 
for N < M — 1, we have 


max ••• max min • • • min v 2 , . . . , 

vj+ieQj+i VNefltf vtelli v j e &j 

(11) < fi -1 (uf,uf;uf,uf;...;u^,u^) 

< min • • * min max ••• max H l {v i, V 2 , . . . , vjv) 

v x eQi Vj€&j Vj+xeQj+i v&eO N 

where 

X i = 1 ? 1 = 1 i • * * ? j 

Xi “1? i — j ~ I* 1 > • • • ? N 

Next we have, N = M and for any t/iefti: 


(^) Xi[H (^i" , it} ; U2 , u 2 ; . . . ; ) 

- ^' 1 (u 1 ,u 1 ;u+,uf;...;uJ f ,u^)] 

< 0 


using the same argument as in (9). 
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Now, for any fixed t>i, t>i; u 2 > • • • ! «m> u m ) is Godunov’s Hamiltonian 

when the initial data for (H-J) 1 has a constant derivative, 

d<po 


dx\ 


(x) = Ui 


Then it follows from (P4) that 


s( _) = T-’. + 5 ( - - 


t ’ 


(where g also depends on vx ). 


By the induction hypothesis, this means we have 

(13) 

< xi&iviM “2 > 

= -Xi5(0,0,... f 0) 

= X 1 X 2 min X 2 " ’Xm min xm H l {v x,v 2 ...v M ) 

view? 

= xi IT 1 (®i,i>2 ,...,um) 

where the extrema is taken on at 02, which depends on ui . The vector 
(t>i, 02, • • • , vn)sQ , i where uiefli is arbitrary. We next take min of the expression 
in (13); If all the x» = 1 or all the Xi = ~ 1 we have equality by (P2). Otherwise 
Xi El, 1 < i < j, Xi — —1? j + l<i<M and we have the right hand inequality 
in (11). Next we have, for any Vj+ifftj+i, following the argument above: 

(14) 

H\u+, Uj ; v % . «2 ; ■ ■ • ; u m- u m) 

> H 1 (u+,u-;...]v j+ 1 ,v j+ x;...;u'l { ,ul f ) 

> Xj +2 min Xj+ 2 ---XM min XAfXi min Xi • • • Xj min Xj ^(vi, *> 2 , • • -»m) 

Vj+ieQj+2 v M eQ M Vjcilj 

We next take the max of the expression in (P4) which gives us the left hand in- 
Vj+l 

equality in (11). 
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We have now obtained formula (11) for any N ; using (4) and (5) give us 
Theorem 1. 

We note that (11) validates the conjecture about Godunov’s Hamiltonian in [7] 
when the inequalities in (13) and 14) become equalities. That paper also discusses 
the high-order accurate non-oscillatory numerical solution of (H-J) in some detail. 

If we take the space gradient of (H-J) and call uj = etc., we 

arrive at the system of conservation laws 

(14) (“*)t + = 0, i = l,...,JV 

o Xi 

with intial data: 

Uj(x,0) = uf if Xi > 0 
= u~ if X{ < 0 
i = 1,... ,n 

Then taking the space gradient of (3) gives us information about the solution to 
this special Riemann problem for a special system of conservation laws. 

We finally remark that if the initial data is convex (concave) or if H(u \ , . . . , u at) 
is convex (concave) then the Hopf formulas [1] for this problem apply. In the case 
N — 1 it was shown in [1] that these formulas give the solution (3) originally derived 
in [6]. (In the one dimensional case the Riemann initial data is automatically convex 
or concave). The same must be true in the multi-dimensional convex or concave 
case. Our general nonconvex results presumably follow from the rather complicated 
formulas in [4], although the connection seems to be unclear. P. Sougandis has 
verified for us (private communication) that the left and right sides of (3) are always 
viscosity sub and super solutions for (H-J). 
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